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Abstract 

We present a simple model where photons propagating in extra-galactic magnetic 
fields can oscillate into very light axions. The oscillations may convert some of the pho- 
tons departing a distant supernova into axions, making the supernova appear dimmer 
and hence more distant than it really is. Averaging over different configurations of the 
magnetic field we find that the dimming saturates at about 1/3 of the light from the su- 
pernovae at very large redshifts. This results in a luminosity-distance vs. redshift curve 
almost indistinguishable from that produced by the accelerating Universe, if the axion 
mass and coupling scale are m ~ 10~^^ eV, M ~ 4 • 10^^ GeV. This phenomenon may 
be an alternative to the accelerating Universe for explaining supernova observations. 
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Current observations of supernovae (SNe) at redshifts 0.3 ^ z ^ 1.7 reveal that they are 
fainter than expected from the luminosity-redshift relationship in a decelerating Universe [||. 
On the other hand, the large scale structure and CMBR observations suggest that the Uni- 
verse is spatially flat, with the matter density about 30% of the critical density It is 
therefore usually inferred that the Universe must have become dominated by a dark energy 
component, which comprises about 70% of the critical energy density, and has the equation 
of state pI p —2/3, implying that our Universe would be accelerating at present. The 
dark energy component could either be a small cosmological constant or a time-dependent 
quintessence energy 0. Neither possibility is elegant from the current vantage point of fun- 
damental theory because of unnaturally small numbers needed to fit the data: the present 
value of the energy density, pc ~ 10^^^ eV^, and, in the case of quintessence, the tiny mass 
smaller than the current Hubble parameter itlq < Hq 10~^^ eV and sub-gravitational cou- 
plings to visible matter to satisfy fifth- force constraints [§. Further, describing an eternally 
accelerating Universe with future horizons is at present viewed as a conceptual challenge for 
string theory and more generally any theory of quantum gravity p|-[TI|. 

Because the SNe observations probe length scales / ~ Hq^ ~ few x 10^ Mpc which are 
inaccessible to any particle physics experiments, it is natural to consider alternative expla- 
nations to the supernova data without cosmological dark energy. A simple such alternative 
is that light emitted by a distant supernova encounters an obstacle en route to us and gets 
partially absorbed |]12|.0 However any mechanism must be very achromatic because the light 
from the SNe appears to be dimmed independently of frequency. This would seem to rule 
out a medium of matter particles, which can absorb light in the optical spectrum but will 
re-emit it in the IR, affecting the CMBR in unacceptable ways. 

In this paper we consider a model where the dimming of SNe is based on flavor oscillations. 
Flavor oscillations occur whenever there are several degrees of freedom whose interaction 
eigenstates do not coincide with the propagation eigenstates. Such particles can turn into 
other particles simply by evolution and evade detection. We will consider a model with 
an axion with a mass m ~ 10~^^ eV, much smaller than the usual QCD axion mass 
scale, 10"^ eV ^ tt^'qcd ~ ^V, but exponentially larger than the quintessence mass 
ttlq ^ 10"^^ eV. This axion couples to electromagnetism through the usual term {a/4M)FF, 
which leads to energy-dependent mixing of the photon and the axion in the presence of an 
external magnetic field B [|T^. Hence light traveling in inter-galactic magnetic fields can in 
part turn into axions, and evade detection on Earth. A source would then appear fainter even 
if the Universe is not accelerating. To satisfy other cosmological constraints we assume that 
the Universe is presently dominated by some form of uniform dark energy which does not 
clump, but need not lead to cosmological acceleration, e.g. with equation of state p/p = —1/3. 

We find that contrary to the familiar example of neutrino oscillations, in our model both 
the flavor mixing and the oscillation length of photons in the optical range are insensitive 
to energy, and so our axion will induce strongly achromatic oscillations of optical photons. 
On the other hand the small axion mass m ~ 10~^^eV insures that the photon-axion oscil- 

* Other possibilities are that SNe may evolve with time or that there are more than four space-time 
dimensions [p^ . 
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lations leave CMBR essentially unaffected. Our task here is to find a small correction to 
the luminosity-redshift relation induced by the oscillations. Because the SNe observations at 
redshifts z ~ 0.5 can be explained by a 10% — 15% total increase in the distance relative to a 
matter dominated universe, the total attenuation of SNe light should be about 30% between 
the cosmic string-dominated geometry and flavor oscillations. Hence we need to account for 
about a 20% of decrease in the luminosity by the flavor oscillations. 
The axion-photon coupling is 

l^int = ^E-B, (1) 

where the scale M characterizes the strength of the axion-photon interactions. This induces 
a mixing between the photon and the axion |[T^,[T^ in the presence of a background magnetic 
field B (as exists in our Universe |18|). Indeed, working in the Coulomb gauge at distances 
short compared to the size of a coherent magnetic domain Ldomi we see that the photon 
with electric field orthogonal to B remains unaffected by mixing. The polarization whose 
electric field is parallel to B mixes with the axion. The field equations are, after rotating 
the coordinate axes such that the propagation is along the ^/-direction. 



where we Fourier-transformed the fields to the energy picture £ and introduced the state- 
vectors I7) and \a) for the photon and the axion. Here B = {e ■ B) \B\ is the averaged 
projection of the extra-galactic magnetic field on the photon polarization e. We will assume 
that the averaged value of B is close to its observed upper limit, and take for the magnetic 
field amplitude \B\ ~ few ■ 10^^ G ||TB,T^|. Therefore its energy density is B^ ~ cH^Mpi, 
where c ~ few ■ 10~^^ and the Hubble parameter is Hq ~ 10~^^ eV. The magnetic fields 
we will be considering are sufficiently small that we can safely ignore the Euler-Heisenberg 
effect |6|,gg]. 



We can now define the propagation eigenstates by diagonalizing the mixing matrix in 
Eq. (^, which is, using B/M = /i. 

This matrix is the analogue of the usual see-saw matrix for neutrinos, with the only difference 
that the off-diagonal terms are imaginary and complex-conjugates of each other. This is 
because the mixing arises from the derivative terms in the field equations rather than the 
potential terms. Defining the propagation eigenstates |A_) and |A+) which diagonalize the 

matrix (^, whose eigenvalues are A=p = ^ =F + A*^^^? we can solve the Schrodinger 
equation (H). The solutions describing particles emitted by a supernova at a distance yo > 
and propagating towards us at ?/ = are 

^ /^^ _ \X_\ g-i[£t+pi{y-yo)] ^f^^ |X \ g-i[£t+p2{y-yo)] 



\a) = , |A_) e-^'^*+P^^^-^°^' + , ^+ |A+) e'^t^^+P^^"^")] , (4) 
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where pk = — \k- It is now clear that as the photon propagates, it mixes with the axion 
by an amount depending on the energy of the particle. In the limit ^ Xi > m? , which 
covers all of the applications of interest to us, the mixing angle is 

sin 6 = I = , (5) 



+ ^ , 
-Ay 



the photon survival probability P^^^ = | (7(^/0) |7(z/))P is 

(6) 



P-y^^ =1 : -^r^ sm 



and the oscillation length is 



AttS 



(7) 



In the limit £ ^ m^/ yU, the mixing is maximal, while the oscillation length is completely in- 
dependent of the photon energy: sin^^ ~ 4=, Lq ~ — • Thus high-energy photons (including 
optical frequencies £ ~ 10 eV as we will see) oscillate achromatically. 

On the other hand, in the low energy limit £ -C m^//i, the mixing is small, and the 
oscillation length is sensitive to energy: sin^ ~ Lq ~ ^f- The oscillations are very 
dispersive, due to the energy-dependence of both the mixing angle and the oscillation length. 
But the probability to find axions P-^^a = 1 ~ ^7^7 is small, bounded from above by 
P^^a < sin2(2^) < V^Vm^ 

In our Universe the magnetic field is not uniform. Assuming that a typical domain size for 
the extra-galactic magnetic field is Ldom ~ Mpc [0,|19|], it is straightforward to numerically 
solve for the quantum mechanical evolution of unpolarized light in such magnetic domains 
with uncorrelated field directions. An analytic calculation shows that in the case of maximal 
mixing, with cos{fiLiiom) > the survival probability is monotonically decreasing: 
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^7-7 - 3 + 3^ 

where the inverse decay length is given by 



^ decay 



In 

For /iLdom ^ 1 this reduces to 

-^decay ~ 



l + 3cOs(^Ldom) 



Thus we see that with a random magnetic field the problem becomes essentially classical 
and after the traversal of many magnetic domains the system is equilibrated between the 
two photon polarizations and the axion. This leads to the generic prediction that on average 
one-third of all photons converts to axions after large traversed distances. 



3 



We can now estimate the axion mass and coupling needed to reproduce SN observations. 
To take the oscillations into account, in the luminosity-distance v.s. redshift formula we 
should replace the absolute luminosity L by an effective one: 

'^e// ~ ^ P-i^^i ■ (11) 

The optical photons must oscillate independently of their frequency. For them, the oscilla- 
tions should reduce the flux of incoming photons by about 20% for SNe at 2; ~ 0.5. This 
requires Ldec ~ /2. Thus the mass scale M for this should be M ~ 4 ■ 10^^ GeV. Note, 
that this is above the experimental exclusion limit for M. The experimental bound on M 
quoted by the PDG is M > 1.7 ■ 10^° GeV ||22|. However, for ultralight axions there 
is p^,|23| a more stringent (though also more model dependent) limit from SN1987A given 
by M > 10^^ GeV, which is still lower than the value required here. 

If the microwave photons had fluctuated a lot in the extra-galactic magnetic field, their 
anisotropy would be very large due to the variations in the magnetic field. To avoid affecting 
the small primordial CMBR anisotropy, AT/T ~ 10~^, the axion mass should be large 
enough for the mixing between microwave photons and the axion to be small. In this limit, 
we can ignore the averaging over many random magnetic domains and simply treat each 
domain as a source of CMBR fluctuation. The photon-axion mixing and the oscillation 
length in that case are given by sin 6^ ~ Lq ~ ^§-- The disturbances of CMBR are 
controlled by the transition probability into axions P-y^a < 4 , which using the explicit 
expression for B is 

^7^^<4-10-^^%P^- (12) 

For microwave photons £ ~ 10^^ eV, and so P-y_»a ^ 2.5 ■ 10^™(eV)'^/m^. Therefore for 
m ~ few X 10"^^ eV we find P^^a ^ 10^'', which is smaller than the observed temperature 
anisotropy. For this mass scale, the oscillation length of microwave photons is Lq ~ \Q~'^Hq^, 
which is of order of the coherence length of magnetic domains Ldom and so a lot shorter than 
the horizon size. This is harmless since the oscillation amplitude is so small. Thus we see 
that if the axion scales are 

m~10^%V, M~4-10"GeV, (13) 

the mixing could produce the desired effect of reducing the flux of light from SNe while 
leaving the primordial CMBR anisotropy unaffected. We stress here that while at early 
times the CMBR photons were much more energetic there were no sizeable extra-galactic 



magnetic fields yet, since their origin is likely tied to structure formation [^. Hence we can 
get a rough estimate of the influence of our effect on CMBR using their current energy scale. 
Having determined the scales, we can check that the approximations we have been using are 
appropriate for optical and microwave photons, respectively. In the former case, the mixing 

4 c 

angle and the oscillation length receive energy- dependent corrections oc ~ 10~ , 

^ optical 

while in the latter case the energy-dependent corrections are oc ^ ^4^^ ~ 10^^, confirming 
the validity of our approximations. 
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Figure 1: The luminosity-distance vs. redshift curve for several models, relative to the curve with 
^tot = (dotted horizontal line). The dashed curve is a best fit to the supernova data assuming the 
Universe is accelerating {^m = 0.3, = 0.7); the solid line is the oscillation model with = 0.3, 
ns = 0.7, M = 4 • 10^1 GeV, m = IQ-^^ eV; the dot-dashed line is = 0.3, = 0.7 with no 
oscillations, and the dot-dot-dashed line is for = 1 again with no oscillations. 



To compare our model with observations, we assume that the constraint on the total 
energy density of the Universe Qtot — 1 is satisfied because the Universe contains some form 
of dark energy which does not clump, but it need not lead to cosmological acceleration. 
A simple example is dark energy with the equation of state w = p/p = —1/3 and energy 
density Qs = 0.7, which could originate from a network of frustrated strings with small 
mass per unit length. Note that because the scale factor a of the universe obeys a/a = 
-j^iptotai + ^Ptotai), and assuming Q^n = 0.3 and Qdark = 0.7, then as long as the ratio 



w = p/p is greater than —1/2.1 ~ —0.48 the Universe would presently not be accelerating. 
These forms of dark energy do not appear to be excluded either by the position of the 
first acoustic peak in the CMBR measurements [^] or by combined CMBR+large scale 
structure fits . In Fig. |I| we have plotted the typical prediction of the oscillation model 
in a spatially fiat Universe with Q^n = 0.3 and Qs = 0.7 against the best fit model for 
the accelerating Universe with a cosmological constant {Qm = 0.3 and = 0.7). The 
two curves are practically indistinguishable. We note that the oscillation model predicts 
limited attenuation of the SN luminosities, unlike some other alternatives to the accelerating 
Universe. The total attenuation is limited to about 1/3 of the initial luminosity, as we 
have explained above. Since for larger values of z the Universe becomes matter dominated, 
and the disappearance of photons is saturated in the oscillation model, the two curves will 
continue lying on top of each other for higher values of z. Thus simply finding higher z 
supernovae [^| will not distinguish between the two models. The main difference between 
the two is that the curve for the oscillation model is an averaged curve, with relatively large 
standard deviations. Therefore it may be much easier to explain outlying events than in the 
case of the accelerating Universe. 
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Figure 2: The effect of the osciUation on the polarization of the photons. We assume that the 
photon emitted is totally polarized, and show how much polarization remains as a function of 
physical distance traveled. 



Let us now consider photons which may pass through the magnetic field of a galaxy, or 
just skim it. The galactic magnetic fields are much stronger than the extra-galactic ones, 
Bg ~ yuG ~ lO'^B. However, the density of baryons (and therefore also of electrons) is 
large enough in such regions that refraction has to be taken into account, which introduces a 
diagonal element Aiu for the photon in A simple estimate |jT6| for this term gives A^n ~ 



10~^^(eV)^ for 10 eV photons traveling within a galaxy, while the off-diagonal terms are of the 
order 10~^^(eV)^. Therefore this term will dominate the mixing matrix, and the oscillations 
will be highly suppressed while passing through the magnetic field of a galaxy. However, 
since there is no evidence for the presence of gas uniformly distributed between clusters, this 
effect is likely negligible for most of extra-galactic space. This is because a simple estimate 
shows, that even assuming the worst-case scenario where all matter is uniformly distributed 
and totally ionized, A^n would be ~ 10~^^(eV)^ for flbaryon = 0.05, which would somewhat 
suppress the mixing. However since matter is not uniformly distributed and definitely not 
ionized in the inter-cluster voids (which make up most of space) this effect should be negligible 
for our results. 

While it is natural to wonder if there are laboratory constraints on our mechanism, a 
simple order of magnitude estimate shows that it would be difficult to observe in a lab. Since 
we have assumed that the extra-galactic magnetic field is ~ 10~^^ T, for a uniform magnetic 
field about 10^^ times larger the corresponding oscillation length would be Lo/10^^ ~ 6- 10^^ 
cm, which is about a thousand times the circumference of the Earth. The current direct 
experimental bounds quoted by the PDG on the coupling of an axion-like particle 



(with a mass less than 0.03 eV) to E ■ B is M > 1.6 x 10^ GeV. 

Another question is whether the oscillations may cause any observable polarization effects 
on the light arriving from the SNe. If the orientation of the extra-galactic magnetic field 



^We thank Georg Raffelt for pointing out this effect. 
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were constant, and the field perfectly homogeneous, light from the SNe would be partially 
plane-polarized. However, since the coherence length of the extra-galactic magnetic field is 
of order ~ Mpc, the direction of the magnetic field is effectively random, and thus no strong 
polarization effects are expected for faraway SNe. Rather, the converse effect of depolarizing 
incoming light is more important, since the oscillations in a random magnetic field may 
deplete existing photon beam polarizations. Because there are distant sources which are 
partially polarized, with the polarization direction correlated with the shape of the source, 
it is important to show that the photon-axion mixing does not completely depolarize light 
from a polarized source. The observed polarization as a function of distance is shown in 
Fig. H, where we see that the polarization decrease is rather slow. This should be expected 
because the degradation of polarization occurs after an axion produced by a polarized photon 
conversion regenerates a photon of a different polarization, after the orientation of the B 
has changed. This is a second-order effect, and so polarization is depleted more slowly than 
intensity. As a result the existing measurements of polarized optical photons from distant 
sources can be accommodated in this model. 

An axion with the scales which are needed for our model can for example arise from 
the spontaneous breaking of an axial lepton number symmetry. Suppose that it couples 
to the electroweak gauge theory in the standard way. Specifically it would couple to the 



electromagnetic field like the QCD axion [15|. Then the mass scale M would be related to 



the scale of the spontaneous breaking of axial symmetry by 

M=-^, (14) 

where a = g'^/4:7T is related to the gauge coupling constant, and ^ is a dimensionless number 
depending on the precise couplings to fermions. We will take a ~ 1/30 and ^ = 0{1) in 
what follows for simplicity. Hence fa ~ 10~^M. In perturbation theory the shift symmetry 
a —>■ a + c protects the axion from acquiring a mass term (more generally any potential). 
This symmetry is broken by nonperturbative effects induced by instantons, which give rise 
to the axion potential [0. Because by assumption our axion couples to electroweak gauge 
fields, a possibility to generate the potential is via the electroweak instantons. In particular 
the axion potential will be of the form 

V{a)=A^[l + cos{^)] . (15) 

Ja 



For example, in a particular supersymmetric (SUSY) model [^], the scale A is 

A^ = e'^^^^e^^MlusyMpi , (16) 

where Msusy is the soft SUSY-breaking mass scale, e a flavor symmetry breaking parameter 
and 02 (Mp;) the electroweak gauge coupling strength at the Planck scale. It is straightfor- 
ward to verify that for Msusy ~ few TeV, e = 0{\) and a2{Mpi) = 1/23, we get A ~ 10 
eV. Since the axion mass is 

A2 

m ~ — , (17) 

Ja 
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we see that for /a ~ 4 ■ 10^ GeV we find m ~ 10^^^ eV. As we have seen above, these are 
roughly the scales most interesting for photon-axion oscillations in extra-galactic magnetic 
fields. 

It is important to stress that while our axion must be light it is not light enough to be 
quintessence. Cosmologically, the axion particles with mass m ~ 10~^^ eV are relativistic 
throughout the history of the Universe, and so would behave like warm dark matter. Because 
they are weakly coupled, with ~ 10^^^ GeV~^, they are out of equilibrium from a very 
early time. If they are not significantly produced during reheating after inflation, their 
abundance can be harmlessly small. On the other hand, the homogeneous axion background 
field a{t) will oscillate around its minimum, with its energy density scaling as cold dark 
matter at late times. Thus one may worry about the cosmological moduli problem which 
such fields usually lead to. However in our case this does not happen because fa ~ 10^ 
GeV and m ~ 10~^^ eV. In the early Universe, the background field will satisfy the slow 
roll conditions, and remain frozen until the Hubble scale comes down to if ~ 10~^^ eV, 
when the Universe cools to the temperature Tj ~ 100 keV. At that moment, the field may 
start rolling. Its kinetic and potential energy will rapidly virialize, after which the energy 
density stored in it will scale as p ~ pi{T /TiY . The initial energy density is determined by 
the initial displacement of the axion from its minimum, which is of order fa- Therefore the 
energy density will be of order given in Eq. ([T6|), A^ ~ (few x 1 eV)^. This would not compete 
with radiation until the temperature comes down to T ~ Pi/Tf ~ 10~^^ eV, which means 
that even if the axion was displaced from the minimum it would remain tiny for a long time 
into the future. Furthermore, while an axionic sector can give rise to both domain walls and 
cosmic strings in the early Universe, because the axion scales in the model we discuss are so 



low, these defects may remain negligible well into the future of our Universe ||30[ . 

In summary, we have presented an alternative explanation of the observed dimming of 
SNe at large distances. The effect is based on a quantum mechanical oscillation between 
the photon field and a hypothetical axion field in the presence of extra-galactic magnetic 
fields. This would result on average in about a third of the photons emitted by distant SNe 
oscillating into axions. This is, roughly, the right amount needed to explain the supernova 
observations. If the average magnetic field is of the order 10~^ G, and the average domain 
size is of order ~ Mpc, one would need an axion whose coupling to the photon is given 
by M ~ 4 • 10^^ GeV, and mass m ~ 10~^^ eV. With these parameters the luminosity- 
distance vs. redshift curve is almost indistinguishable from the curve of an accelerating 
Universe with fi^ = 0.3, f^A = 0.7. Since the precise value of the luminosity-distance for a 
particular supernova depends on the details of the inter-galactic magnetic field, we expect 
more variations in the observed luminosity, and thus this model may easily incorporate 
outlying data points. However, distinguishing this model from the accelerating Universe 
paradigm will likely be easier through improving the bounds on the couplings of ultra-light 
axions, by understanding the details of the intergalactic magnetic field, or by a precise 
independent determination of the equation of state for the dark energy component, for 
example through the DEEP survey pT|. 
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